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O- ■ Abstract 

o 

Within the Closed Time Path (CTP) framework, we derive kinetic equations for 
particle distribution functions that describe leptogenesis in the presence of several 
lepton flavours. These flavours have different Standard-Model Yukawa couplings, 
which induce flavour-sensitive scattering processes and thermal dispersion rela- 
tions. Kinetic equilibrium, which is rapidly established and maintained via gauge 
interactions, allows to simplify these equations to kinetic equations for the ma- 
trix of lepton charge densities. In performing this simplification, we notice that 
the rapid flavour-blind gauge interactions damp the flavour oscillations of the lep- 
tons. Leptogenesis turns out to be in the parametric regime where the flavour 
oscillations are overdamped and flavour decoherence is mainly induced by flavour 
sensitive scatterings. We solve the kinetic equations for the lepton number densi- 
ties numerically and show that they interpolate between the unflavoured and the 
fully flavoured regimes within the intermediate parametric region, where neither 
of these limits is applicable. 



1 Introduction 

Calculations of the baryon asymmetry of the Universe require the description of CP- 
violating processes in a finite- density background. Conventionally, semi-classical Boltz- 
mann equations are employed, where one computes the evolution of classical particle 
number densities in the presence of scattering processes, that are given by matrix el- 
ements from the in-out formalism. Since CP violation in these matrix elements is a 
higher order effect that is typically induced by loop diagrams, where particles in the 
loop are kinematically allowed to be on shell, this approach bears notorious difficulties 
concerning the correct counting: Is a certain process higher order or is it a combination 
of two leading order processes, that are already accounted for in the Boltzmann equa- 
tions? For the CP-violating processes in baryogenesis, a way of addressing this caveat 
is the method of real intermediate state (RIS) subtraction [lj. 

A more direct and systematic way to deal with this counting problem is to avoid the 
detour via the in-out matrix elements and instead to formulate the problem directly in 
terms of Green functions. A method that achieves this is given in terms of the Closed 
Time Path (CTP) or in-in formalism |2|,|3], which has been successfully applied to various 



models of baryogenesis including also a realistic scenario of leptogenesis [4j-|l0| . The set 
of Green functions that are computed in the CTP framework encompasses both, the 
spectral and the statistical information of the system. To leading accuracy in a weakly 
interacting situation, it is often sufficient to approximate the spectral functions by on- 
shell (^-functions, while the statistical contributions encode the distribution functions of 
quasi-particles. In this way, the kinetic equations that govern leptogenesis have been 
recovered systematically and additional corrections due to the full quantum statistical 
distributions of right-handed neutrinos and the leptons and Higgs boson of the Standard 
Model have been derived J9|. 

The discussion in Ref. [9| is for a single flavour of Standard Model leptons £, which is 
appropriate in those situations where the different Standard Model Yukawa couplings of 
these flavours are negligible, because an asymmetry is only produced for one particular 
linear combination of the lepton flavours. However, at smaller temperatures, when the 
lepton Yukawa couplings h of the Standard Model approach equilibrium, the flavour 



degeneracy is broken and effects of flavour have to be taken into account |llHl4| . The 
derivation of a set of kinetic equations from non-equilibrium quantum field theory, which 
covers both the unflavoured and fully flavoured regimes, and is valid in between, is still 
missing up to now, and is the subject of this paper. 

A simple model that encompasses the salient features of flavoured leptogenesis and 
which we consider in the present work is specified by the Lagrangian 

C =-^ N i(i$> - Mi)ij Nt + i> ia i$i) ia + t/W#r 6 + (^0 t )(9 M 0) (1) 
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where ipNi is the four-component Majorana spinor representing the right handed singlet 
neutrino N t , ip^ a is the spinor for the SU(2)l doublet of left handed Standard Model 



leptons £ a , where a is the flavour index, and ip^ a are the corresponding charged right 
handed leptons. The SU(2) L doublet of Higgs fields is represented by 0, and we define 4> = 
(e</>)T with e the antisymmetric 2 matrix in the SU(2)l indices with e 12 = 1. For simplicity, 
we assume i = 1,2. The generalisation to more than two right handed neutrinos iVj is 
straight forward . 

Single-flavour calculations for leptogenesis can be easily generalised to the multi- 
flavour situation, provided there is a flavour basis in which one can express the lepton- 
number densities of the several flavours and at the same time neglect possible correlations 
between the different flavours. When the lepton Yukawa couplings are fully in equilib- 
rium, the appropriate flavour basis is where these couplings h are diagonal. In contrast, 
in the unflavoured regime the appropriate basis is determined by the linear combination, 
in which the lepton asymmetry is produced. 

These considerations raise the following questions: First, is there a kinetic equation 
that is manifestly covariant under the choice of the flavour basis? And second, can 
such an equation also deal with the intermediate regime where the couplings h are not 
yet in full equilibrium but are non-negligible at the same time? The latter point is of 
particular importance since leptogenesis is a process that takes a finite amount of time 
while temperature is decreasing, such that the flavour-sensitive interactions may be out 
of equilibrium initially, but equilibrate at later times. Clearly, in order to address this 
question, one has to promote the set of lepton number distributions within a certain 
flavour basis to a matrix, that also allows for off-diagonal flavour correlations. In a 
heuristic approach, by appealing to the Hamiltonian evolution of a density matrix, such 



a set of equations has been proposed 13(, and numerical solutions to these equations 



have been obtained 15 



On the other hand, it is clear that the Green functions within the CTP formalism 
allow for the possibility of off-diagonal correlations in a straightforward way. Indeed, once 
the appropriate Kadanoff-Baym equations for unflavoured leptogenesis are known [9|, 
the multi-flavour generalisation is easily written down. However, the Kadanoff-Baym 
equations are a coupled set of integro-differential equations that needs to be subjected 
to some analytic approximations before solving it numerically. Key simplifications are 
first the gradient expansion, which effectively is an expansion in terms of time derivatives, 
in deviations of the distribution functions from equilibrium, and in the small coupling 
constants. Second, one makes use of a separation of time scales between the interactions 
induced by the small Yukawa couplings Y and h and the faster gauge interactions, that 
maintain kinetic equilibrium and induce lepton-antilepton pair creation and annihilation 
processes. This allows to describe the distribution functions by generalised chemical 
potentials. 

The application of these strategies of simplification to the multi-flavour case con- 
stitutes the main body of the present work. In Section [2j we set up the multi-flavour 
Kadanoff-Baym equations and show that the flavour-dependent dispersion relations for 
the leptons, that are induced by the right handed neutrino and Yukawa couplings Y 
and h, respectively, give rise to commutator terms in the kinetic equations which are 
also characteristic of flavour oscillations in the presence of tree-level mass terms. In 



Section |3l we make use of the short time scale for kinetic equilibration in order to ex- 
press the lepton densities in terms of a generalised chemical potential. We show that the 
gauge interactions in addition to enforcing kinetic equilibrium damp the flavour oscilla- 
tions and that parametrically, flavoured leptogenesis is in the overdamped regime. The 
processes that turn out to dominate the decay of the off-diagonal flavour correlations are 
the flavour- sensitive three-body scatterings through the couplings h. 

Let qiab denote the Hermitian matrix of lepton charge densities, including flavour 
off-diagonal correlations, which is defined more precisely below. The main result of the 
present work is the kinetic evolution equation 

—^ — / , [Qiac^cb — ^aclicb ~ W ac q£ c b ~ q.(,a(Wcb\ + 2S a b — T £ab . (2) 

dr) V 

In this equation, rj denotes the conformal time, which is related to the physical comoving 
time through dt = a(rj)dri, where a(rj) is the scale factor of the expanding Universe. The 
matrices W, S and Tf are all Hermitian. Lepton number and CP-violating source terms 
are represented by S, W encompasses the washout rates for the various lepton flavours, 
and Tf is the matrix of flavour-sensitive damping rates. In case we impose that qi a b is 
evaluated in the basis of mass eigenstates of leptonic quasi-particles, the anti-Hermitian 
matrix S compensates for possible time-dependent flavour rotations. Numerical solutions 
to this equation are presented in Section HJ and it is shown that it indeed interpolates 
between the fully flavoured and the unflavoured regimes. 

We conclude in Section |5j Appendix |A] contains a discussion of the constraint and 
mass-shell equations including finite width effects, and the contribution of the right- 
handed neutrino to the thermal mass of the leptons is calculated in Appendix [B] 

2 Flavoured Leptons 

The usual strategy for deriving kinetic equations is to decompose the Schwinger- Dyson 
equations on the CTP into equations for the retarded and advanced propagators and 
the Kadanoff-Baym equations. In the weak coupling limit, one can solve the latter when 
approximating the particle densities as purely on-shell. This corresponds to taking the 
particle width to zero. 

In this Section, we reiterate these arguments for the left-handed leptons. As an 



extension of earlier approaches |16Hl8l|. where mass terms are introduced at tree- level, 
we describe here the dynamics that arises from thermal one-loop corrections to the 
dispersion relations, which are mediated by flavour-blind gauge interactions as well as 
flavour-sensitive Yukawa interactions. 

2.1 Schwinger-Dyson Equations 

We employ here the notations and conventions for the CTP formalism and the gradient 
expansion that are explained in more detail within Ref. [9|, [l9|, |20j • The Schwinger-Dyson 



equation for the flavoured left-handed lepton propagator is 

iMfi(«, «) = /*''<W*(« - t;)P R + J2 I d^tad^ ™)<*K «) , (3) 

ft ^ 

which can be more compactly written as 19|, |20| 

ys{" = fSf'SPn + J2$ h ® s i 9 > ( 4 ) 

where the symbol denotes a convolution and <5 the delta function, f,g,h = ± are 
the CTP indices and a, b, c flavour indices. The CTP structure can be decomposed into 
equations for the retarded and advanced propagators and a Kadanoff-Baym equation, 

[@S?' R = SP R + Sf Sf' R ± i# Sf' R , (5a) 

i^<'> = Jf © S<-> + ^ Sf + \ (E> 5? - S £ < © Sf ) . (5b) 

Here and in what follows, we make use of the definitions 

G< = G + - , G> = G- + , G T = C++ , G f = G- (6) 

and the combinations 

G A = G T -G> = G< - G f (advanced) , (7) 

G R = G T -G < = G > -G f (retarded) , 

G H = -(G R + G A ) (Hermitian) , 

G A = -(G A - G R ) = -{G > - G<) (anti-Hermitian, spectral) , 

where G may stand for any two-point function on the CTP, in particular the propagators 
and self-energies. 

The Wigner transformation is defined as a Fourier transformation of a two-point 
function with respect to the relative coordinate, while keeping the average coordinate 
fixed: 

d 4 r e lk ' r G(u, v) , where r = u — v and x = — - — . (8) 

It provides a separation between the typical energy scale (which is the temperature T in 
the present case) and the small macroscopic inverse time-scale, that governs the evolution 
of state parameters (particle number and charge distributions in the present case). The 



ratio of these two time scales is then used to define the gradient expansion. In Wigner 
space, Eqs. (JSJ) take the form 



e" 10 



-tf=Ftf\{st>*\=Pk, (9a) 

{t-tf}{Sf> > }-e^&f>}{Sf} = ±e-»({tf}{Sf} -{$<}{#}) , (9b) 



where 

o{A(k,x)}{B(k,x)} = l - {[d^ x) A(k,x)]d {k)fi B(k,x) - [d {k) ,A(k,x)}d^B(k,x)) . (10) 
Eqs. Qj correspond to an infinite tower of integro-differential equations. Approximate 



solutions can be obtained within the scheme of gradient expansion [21| . In the context of 
leptogenesis, gradients occur due to the deviation of particle number distributions from 
equilibrium, which is induced by the Hubble expansion of the Universe. At the same 
time, for a sizeable lepton asymmetry to occur, it is crucial, that during leptogenesis, the 
rate |Yi a | 2 T is not very different from expansion rate H . Similarly, within this paper we 
are particularly interested in the parametric region where the r-lepton Yukawa coupling 
h T relates to H as ti^T ~ H, since otherwise we are either in the fully flavoured or 
in the unflavoured regime, which can be described by conventional approaches. In the 
present context, gradient expansion is therefore understood not only as an expansion in 
time-derivatives, but also as a perturbative expansion in Y and h. 

Performing the expansion of Eqs. ([9]) up to first order in gradients, we obtain 

Sf' R = PR, (Ha) 



Tt T HJ 



±?S<>> + $- yf)Sf'> - fr>S? = \ $>S< - $fS>) . (lib) 

The derivative term ^$S e ' would contribute to Eq. (II lap only at second order in gra- 
dients, since the retarded and advanced propagators do not depend on the particle dis- 
tribution functions at tree level. Damping occurs explicitly through the collision term 
on the right hand side of Eq. (Illbj) (for a detailed discussion, see Ref. (22|). At the same 
time, damping is contained in the fy term in Eq. (II lap , as it is well known from linear 
response theory. 

In order to take account of the dilution of particles due to the expansion of the 
Universe, we follow Ref. [9[. We first observe that after appropriate field redefinitions 
and the rescaling Mi — > a(rj)Mi, the Lagrangian (JI|) describes the fields in the background 
of a flat Friedmann- Robertson- Walker Universe in conformal coordinates defined by the 
metric 

S> = a 2 (r]) diag(l, -1, -1, -1) . 

We then take all explicit momentum variables within the equations for the Wigner 
transformed quantities as conformal. The relation to a physical momentum is given 



by /c p h = k/a(rf) and the time coordinate is understood as the conformal time rj. Like- 
wise, T denotes a comoving temperature that is related to the physical temperature as 
T p h = T/a(r)), and = 1/T. The masses of the right-handed neutrinos M; are the 
physical masses. When they occur explicitly in the Wigner transformed equations, they 
are accompanied by the scale factor 0(77) to give the conformally rescaled mass a(rj)Mi. 
The Hermitian and anti-Hermitian parts of the Kadanoff-Baym equations ( 111b|) lead 
us to the constraint and the kinetic equations 

2k°i 1 °S<'> - {k ■ 77 ° + ^V,i 7 X'>} - {i^> 7 ,7°Sf } = -\ fa ~ iCj) , (12a) 
id v ij S<'> - [k ■ 77 ° + ^ 7 °, h°S^] - [i^'V, l°Sf] = -\ (iC t + iC|) , (12b) 



with the collision term 



C e = i^iS< - iftiS? . (13) 



2.2 Thermal Self Energies 

We now specify the form of the Hermitian self energy ^ , which determines the thermal 
corrections to the dispersion relation. For the purpose of this discussion, we assume that 
the deviation of the right handed neutrino distribution from thermal equilibrium is small, 
so that the self energy, being proportional to coupling constants, can be evaluated with 
thermal propagators to first order in the gradient expansion. The self energy receives 
contributions from the interactions specified in Eq. (pQ), but also from SU(2)l x U(l)y 
gauge interactions. We parametrise this as 



Pi 



K 



7 V + T fl ) + ^ (?" + «* ~ si Mk°) [? bl + <T fl ]) 



Pl , (14) 



where we have decomposed the contributions to the self energy into a flavour blind part 
that originates from SU(2)l x U(l)y gauge interactions 

^l(k°, k) = S ab ^ l (k°, k) , ^(A; , k) = 6 ab , h \k°, k) , (15) 

and a flavour dependent part that receives contributions from the charged lepton and 
the singlet neutrino Yukawa couplings. To one loop order, these can be parametrised as 

&(*°, k) = hih cb ^ h (k°, k) + J2 y: a y^ Y (k\ k) , (ie) 

i 

£(*°, k) = hlh^'Hk , k) + J2 OW^* , k) . 

i 

In the hierarchical mass limit (Mi <C M 2 ), which we assume within this paper, we may 
restrict the sum to i = 1. Note that 7°^ is Hermitian, such that ? fl and ^ fl are Hermitian 
matrices in flavour space. 

6 



The fact that two independent functions s(k°, k) and <;(k°, k) occur is because the 
self energy ^ acquires contributions that are proportional to ^ and ^, where u^ - 



(1, 0, 0, 0) T is the plasma vector 23|. The relation between s(k°, k), s(k°, k) and the terms 
proportional to ^ and ^ can be easily established, see e.g. Eqs. (IB.2[) . The motivation 
for our parametrisation is that it corresponds to a decomposition into a correction <; for 
the dispersion relation and a correction ^ that leaves the dispersion relation unaltered 
[see Eqs. (1261) and (1A.8I) below]. Besides, <^ and ^ exhibit useful symmetry properties 
under the exchange k° — >■ —k° [see Eq. (|35|) below]. 

The matrix <^ fl can be diagonalised through a unitary transformation U and we define 



*S 



tfW. 



(17) 



Note that in general, [/ is momentum- and time- dependent. At temperatures M 2j 3 ^> 
T/a(rj) ^> Mi and momenta |k| ~ T, both <^ fl ' h and <^ fl ' r are approximately proportional 
to T 2 /|k|, such that the diagonalisation matrix U is constant in time. When the tem- 
perature T/a{rj) drops below Mi, the distribution of N\ becomes Maxwell suppressed. 
We consider this situation in Appendix [Bl The function ^ ' then falls toward zero as 
~ [T/(aMi)] 4 , cf. Eq. (IB. lip , and U may generically undergo a change. Afterwards, 
at temperatures T/a(rj) <C Mi, U becomes constant again and the matrix Wh/hU is 
diagonal. 

All quantities that carry left-handed flavour indices transform under the basis trans- 
formation defined by U. We denote matrices evaluated in the flavour- diagonal basis by 
a subscript D. For example, 



Hd — U Tt u > 



(18a) 
(18b) 



Note that unlike <^, these matrices are in general not diagonal in flavour-space. 

Inserting the definitions (118aj) and (118bj) into (112bl) and multiplying with W from 
the left and with U from the right, we obtain the kinetic equation in the lepton thermal 
mass basiaj 



i«V 7 °^> + i [3, h°S< 6 >] - fc, 7 °, i 1 °S l 



'id 



1 



iC 



ID 



\C 



CD 



where Cm = WCelf and 



S = tfdrjU 



(19) 



(20) 



is the compensation matrix for time-dependent flavour rotations. At first order in gra- 
dients, the only consequence of a time dependent U is the additional term involving S. 
Therefore, we switch to the diagonal basis and drop the subscript D from all subsequent 
expressions. 



1 The terms 77 and 7°S'/' > commute when using the ansatz (J2TJ) below. 



2.3 Kinetic and Constraint Equations 

The Weyl fermion propagator can be parametrised through the vector and pseudovector 
functions 

W* = \ E W i 1 + h * • 757 ° 7 ) + ^ ( 75 + h * • 7 ° 7 )] ' (21) 

h=± 

where k = k/|k|. When compared to the case of a Dirac fermion, there are no scalar 
and pseudoscalar contributions. This is because gauge symmetry prevents the dynamical 
generation of scalar, pseudoscalar and tensor densities, provided the gauge symmetry is 
neither broken spontaneously nor through initial conditions, as we assume here. Thus, 
Eq. (T2T1) is the most general form of the lepton propagator compatible with isotropy 
and chiral symmetry. Besides, from the fact that the leptons £ are left-handed, we 
immediately obtain 

9tt = 9nt - 9^ ■ (22) 

We furthermore see that k ■ 77 and i7°S'/' > commute, such that the constraint and 
kinetic Eqs. (1121) simplify to 

2(k° - k ■ 77 °)i 7 °S<>> - {^ 7 °, i 7 °^>} - {i#- V, j>Sf) = -\ (iCt - iCj) , 

(23a) 

id^S^- [j?/,i7°S<'>] - [i^'V.ySf] =-|(iC, + iC]) . 

(23b) 

To zeroth order, the constraint equation f)12ap reduces to the simple form 

{lf,iS<> > }=0. (24) 

Substituting the ansatz (T2"T1) leads us to 

g< > k° + h\k\g< 3 > = 0, (25a) 

^ < 3 > A; + ^|%< > = 0. (25b) 

The constraint (12^|) then implies that g^' > (k°, k) is non-vanishing only when k° = — h\k\, 
which corresponds to the singular zero-mass shell. In particular h = —sign(k°), that is 
for leptons (A; > 0) the helicity h = — 1 is negative, while for anti-leptons (k° < 0) 
the helicity h — 1 is positive, as expected. This relation is weakly broken for momenta 
|k| ~ T when including thermal corrections, because hole modes exhibit an opposite 
connection between frequency and helicity. For momenta |k| ^C T, the hole modes 
couple to the plasma at a similar strength as the particle modes. However, this region 



only corresponds to a small portion of the available phase space, such that we may 
neglect it here. 

Substituting the parametrisations (I14p . (I2ip and the constraint (122]) into the kinetic 



equations (123b|) and taking the trace leads us to 



id v g^ + i [S, g<>] + h [<r fl , g<'>] = -itr (ifc + iCJ 



(26) 



where the trace is taken only in Dirac space, and ^ fl and all other objects are evaluated 
in the mass-diagonal basis. As explained above, the helicity is determined from the 
zeroth order constraint equation by the relation h = — sign(/c°). Eq. ( 126]) is accurate 
up to first order in gradients, because <^ fl itself is of first order. We note that the term 
[i^ =,> 7°, 7°S'f i ] in Eq. (]23bj) does not contribute to Eq. ([2"6"j) at first order, since first, 
Sf can be evaluated at zeroth order, where is it independent of the particle distribution 
functions and therefore proportional to the unit matrix in flavour space; and second, in 
Dirac space the commutator reads 



[iJJ'V.Yttf] oc-i[P R (a<>> ^ + 6<'> 7 )PL7 ,7°i 3 L^R] =0. 



(27) 



Since we work in the flavour-diagonal basis, where <^ fl is diagonal, the commutator in- 
volving q in Eq. (126 j) can be explicitly evaluated, which yields 



®v9t£ + i [H, 9^] ab + h(C - <&)g 



hab 



-~ ( 

4 



iC f + \C 



<ib 



(28) 



Hence, the thermal dispersion relations have the same im pac t on the equation of motion 
for the lepton density as explicit Dirac masses would 16M18 . 



An important feature of Eq. (126]) is the sign change of the commutator term involving 
the thermal dispersion relation through <^ fl when h — > —h or, alternatively, A; — > —k°. 
We now show that at the one-loop level, the elements of h<^ = — sign(fc°) <; R (k°, k) are 
indeed odd functions in k°. We define charge and parity conjugation through 



rfP{x) = C^p T (x) 
i/j p (x) = Ptp(x) , 



(29) 
(30) 



where u = (u , — u) and where in the Weyl representation, the conjugation matrices are 
given by C = i r j° r j 2 and P = 7 . Thereby, we fix possible CP phases that can arise in the 
definition of these conjugations to zero. The charge and parity conjugate propagators 
are 



U3%'(v,v) =<^2(« / )^S(«')> = C \iS! b {(v,u 

«(n f )^b(v 9 )) 



& 



: p J9< 



iS iab^ U 



:fg 



PiSZ{u,v)P^ 



(31a) 
(31b) 



The transposition acts here only on the Dirac indices, which in contrast to the flavour 
and CTP indices are not written explicitly. The CP conjugate of the Hermitian self 
energy is then given by (we suppress the average coordinate in the argument of q and q) 



P,H 

ab 



(k,x) 



CP 
Pr 



7f ba (-kx)] (cpy 



(32) 



+ riga(*°)[fi(-fc) + £(-*)]) 



Pi. 



On the other hand, we may calculate this self energy from the CP-conjugate Lagrangian, 
within which the coupling constants are complex conjugated, as 



yCP,H 

Hab 



(k,x) 



(lb 



(k,x) 






(33) 



provided the initial conditions preserve CP symmetry (no primordial asymmetry). To 
the one-loop order, the coupling constants appear as the prefactors h'h and Y*Y within 
fa? , cf. Eq. (fT5"j) . The effect of CP conjugation therefore amounts to the replacements 
[tfhU -> [tfh]* ab = [tfh} ba and [yty] a6 _+ [yty]* 6 = [Y^Y] ba , and it follows that to 
one-loop order 



Jtab \k,x) — ff ba {k, 



x 



Comparing this to the relation (|32|) and substituting Eq. (Tl4|) . we find that 



d( 



-k°, k) = -q* b (k°, k) and q! b (-k°, k) = q* ab (k°, k) 



(34) 



(35) 



and accordingly for q and q , which implies the result to be shown for the elements of 
^ fl . These symmetry properties with respect to k° are in accordance with the equilibrium 



results from Ref . 23| and Appendix [Bj The present argument shows moreover that they 
also hold under out-of-equilibrium conditions. 

For calculating the momentum integrals in the collision terms, we use on-shell condi- 
tions that we obtain from the constraint equations. In order to achieve accuracy to first 
order in gradients, it again suffices to solve the constraint equations to zeroth order, since 
the collision term is suppressed by higher orders in the coupling constants. A similar 



approximation is applied in Ref. [18|, where in contrast to the present work, within the 
constraint equations, small tree-level mass differences rather than differences between 
one-loop dispersion relations and finite widths are neglected. A general solution to the 



zeroth-order constraint equation (124)) is given by the Kadanoff-Baym ansatz 
iS& = -2Sf [$(k°)f+ b (k) - $(-k°)(l ab - ffj-k))] , 



iS' > 



-2St [-#(k°)(l ab - /+,(k)) + i)(-k°)fi ab (-k)} 



(36a) 
(36b) 
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where 



Sf = nPaPnS(k 2 ) 



(37) 



and f lah are the distribution function matrices of leptons and anti-leptons. Comparing 
Eqs. (US} to Eq. (J2U) we identify 

g< = -2n6(k?)0(k o )W+(k), (38) 

g> = 2n5(k 2 )#(k°)\k\(l - f+(k)) , 

g< = -2n5(k 2 n-k )\k\(l - fr{-k)) , 

gX = 2<K5(k 2 )$(-k°)\k\ff(-k), 

which is consistent with h = — sign(/c°) and k = ±|k|. While for the present work, these 
zeroth-order solutions to the constraint equations are sufficient, we show in Appendix [A] 
how to extend them to first order in consistency with the equations for the retarded and 
the advanced propagators (II lap . 

Note that if we identify ff ah with expectation values of number density operators 
then it follows from Eqs. (1361) and the operator definition of iSf^ that f^ ab ~ (o&Oa) 
corresponds to the lepton density matrix, while f^ ab ~ {b\b h ) corresponds to the transpose 
of the anti-lepton densities. This may also be seen when relating the CP conjugate to 
the original propagator as 



lb. 



CPJg 



(k,x) 



cPhsHi 



-k,x) 



(CP? 



(39) 



E 

h=± 



9hba{ 



(-k) (l- 7 V-(l-7 5 )^k- 7 



'iSebai 



-k,x) 



where we have used the ansatz (121]) and have assumed spatial homogeneity, g 9 h (k) = 
g 9 h (k). Therefore, a sign flip in k yields the negative of the flavour- and CTP-transposed 
CP conjugate propagator. The use of Sf' > (k,x) with k° < rather than S e ' ,:> (k, x) 
with A; > to describe the anti-lepton densities has the advantage that the resulting 
kinetic equations are flavour covariant. This is because the former propagator has the 



same flavour transformation properties as S^' > (k,x) with k° > 0, while the latter trans- 
forms in the complex conjugate representation. This prop erty has been described and 
used before in the context of electroweak baryogenesis 



16, 17 . 



3 Kinetic Equations for Lepton Number Densities 



In this Section, we perform simplifications of the kinetic equations ([28]) . such that they 
attain a form which can be solved numerically. The key simplification arises from the 
separation of the time scales of kinetic equilibration and flavour-sensitive interactions. 
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(2tt) 3 
d 3 k 



Because the former is much faster, the distribution functions are driven to kinetic equi- 
librium, such that they can be approximated by the Bose-Einstein or Fermi-Dirac form, 
parametrised through a matrix of chemical potentials. An integration over the momen- 
tum then allows to express the kinetic equations in terms of charge densities. We shall 
often use spatial homogeneity to set f^ b {— k) = /^(k). 

3.1 Matrices for Lepton Number Densities 

The lepton number density matrices are defined as 

A ' &(k) = - / jgj [ f^tr [WSZ J , (40a) 

- - j w^ ] = / Sf L ^ tr ™ • (4ob) 

Due to the presence of the fast gauge interactions, we may assume kinetic equilibrium 
for the leptons, and denote 5nf ab = nf ab - n^ b q and SfjL = f± b - f^ q . Then, we can 
introduce a matrix of generalised chemical potentials /j, ab for particles and antiparticles, 
and write 

f^)=( em -l ± + 1 ) ab > (41) 

such that the number density matrices are, to first order in the chemical potentials, 

Ka b = /*£^ ■ (42) 

This allows to relate the lepton number densities to the distribution functions: 

*/i(k) = 12^ ^ . (43) 

We also introduce the deviation of the <, > propagators from equilibrium, 

i5S eab = -2Sf m°)Sf+(k) + $(-k°)5f- b (-k)] . (44) 

Besides, we use corresponding approximations and expressions for the right-handed lep- 
tons of the Standard Model by replacing £ — > R. 

The ansatz ( 14TT) is valid provided the interactions that establish kinetic equilibrium, 
which are the pair creation and annihilation processes and scatterings with gauge bosons, 
are faster than processes that distinguish between the particle flavours, in particular 
flavour oscillations and flavour-sensitive damping rates. Provided these assumptions 
hold, which is verified in Section I3.3[ the time scales for kinetic equilibration and for 
flavour effects separate. Since ^ is Hermitian, it is always possible to bring Eq. (J4"T|) 
to diagonal form by a flavour rotation. Due to the separation of time scales, kinetic 
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equilibrium in this diagonal basis is then attained at a rate that is faster than the 
flavour effects that may change the flavour orientation of fi a b. 

Using the ansatz (1361) and the spectral function (13TJ) together with the decomposi- 
tion f[2~Tj) and the constraint (1221) gives the relations 



00,0 

.0 



y?(k) = T2 / l^^k). (45) 

0,-00 

Performing a fc°-integration of Eq. (|28|) and using Eqs. ( 140]) then yields 

00,0 
1 f dk° 

djLM = [2, tfwu t i(d - 4)&(k) ± 2 tr y ^(^ + C L) • (46) 

0,— 00 

Integrating over three momenta and applying the expansions explained above, we find 

00,0 
1 /* dk f d 3 k 

dM ab = [2 cff , SnfU T iALJ&Sni* ± - tr J -£ J j^{C lab + C\ ab ) , (47) 



0,-00 



where we have used that the equilibrium distributions are diagonal in flavour, n e ° b q = 
Sabn^a- We have defined here the thermally averaged frequencies of flavour oscillations 
and the compensation matrix as 

/d 3 k 19/9 3 p' 8 l k l 
7^(^Tiy ( ^ (|k| ' k ' r/) ~^ (|k| ' k ' r/)) ' (48a) 

„ eff . . f d 3 k 12/3Vl k l „,„ . , , , AQU . 

-^ = J(2^W^W^ klKr]h (48b) 



and used that <j fl and S are symmetric in their first argument, cf. (135 j) . The dominant 
contributions to the phase s pac e integrals originate from regions where |k| ~ T, where 
Sab can be approximated by [23] 

&(*°, k) = %f^ + £ O*,^ , k) . (49) 

i 

While the form of Y' f Yq i ' is more complicated in general, as we discuss in Appendix |B| 
it is of the same order or smaller than Ivhc, ' when flavour effects are important. We 
can therefore estimate 

Auf ab = O (h%T) , (50) 

where h T is the r-lepton Yukawa coupling. 
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We decompose the collision terrro as 

Ci = Cj + Cf + Cf . (51) 

The term Cj describes decays and inverse decays of Ni and hence the washout and the 
CP-asymmetric source of the lepton densities. Flavour sensitive interactions mediated 
by the Standard Model Yukawa couplings are encompassed in Cf, while flavour-blind in- 
teractions mediated by gauge couplings are taken account of within Cf . In the following, 
we show that these particular contributions can be cast into the form 

dSn tab _ W e&! ± _ r ± ^cff • a cff r ± (^ 

Q - -ac 0n icb 0n £ac-cb + 1L1U! £ab 0n eab V 0Z ) 

- Ys^Mcb + Snf; a W b * c } ± S ab - T h \6nj ab + 5nJ ab ) - Y^ b . 

c 

3.2 Source and Washout Term 

The washout term for Sn^ ab is given by 

d 3 k f°° dk n 1 



^Iw^r^l' 1 ^^ms^k) ~ VLWS? cb (k)] 



evaluated to order Y? a (certain higher order terms are accounted for by the source term). 
Close to equilibrium, we can write 

i?L(*)iS&(*) - ^LWiS> b (k) = -i $<JA;) - %IM i6S M (k) , (54) 

and we note that 

(27r) 3 2y / k' 2 + {a^M^ 2 ^ Zn > ^ k I 

(55) 

x sign(fc )iW + a^M^P^ (f m (k f ) + f^k")) . 
Substituting Eqs. (143]) and (jSJ), we can identify 

1 f d 3 k d 3 k' d 3 k" 



W ac =- Y*Y lc I - " " , = ' u = - " '" , (2?r) 4 (? 4 (fc / - k - k") 

2 1C i(2vr)32|k| (27r)3 / k , 2 + (a(r/)Mi)2 ( 27 r)32|kf 



x2fc.fc'(/ Jvl (k') + /^(k")) 



(56) 

12/ 33 e /3|k| 



fe^l + r 



2 The definition of the coliision terms C in Ref. [9( differs from the present ones by an additionai 
integration J dk° /(2-k). 
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The washout term for 5nJ ab follows correspondingly. 

In straightforward generalisation of the single flavour case, the CP-violating source 
term for Sn^ is 

s * = I i E KYZYicY* - Y 2 * a Y 2 * c Y lc Y lb ] (57) 

c 

MA f d 3 k> ^(kQS^kQ 

TT / / 0/ati(kJ, 

M2 / y (2vr) 3 2y / k' 2 + (a^MO 2 ^ 

where l9| 



^) = &» / (27r y 2 | P | (2 ^2| q | w 4 ^ k -p - *)p* ft - /r(p) + /; q (q)) , (58) 

with g w = 2, and where the source for the anti-leptons SnJ ab is —S a b- The deviation 
of the distribution function of the right-handed neutrinos from equilibrium is denoted 
by <$/jvi(k) = /jvi(k) — /^(k). This source term is understood to include both, wave- 
function and vertex contributions in the hierarchical limit, M x <C M 2 . Note that there is 
an additional wave-function contribution that violates lepton flavour but conserves total 
lepton number ll|, ll4|. Within flavoured models of leptogenesis, this may contribute to 
the final lepton asymmetry. Compared to the total lepton-number violating contribu- 
tions, it is however suppressed by a factor Mi/M 2 since one picks up the ft ~ M 1 rather 
than the M 2 term from the numerator of the intermediate neutrino propagator in the 
wave-function diagram. Hence, we do not account for this term here. 

3.3 Flavour Blind Interactions 

The flavour-blind contribution to the lepton self-energy that is mediated by gauge inter- 
actions can be expressed as 

& =g2 1 WYW)i {2nYnk ~ k '- w®&wi^%*m , (59) 

where iA^ is the gauge boson propagator on the CTP. The corresponding collision 
term is then 

C(fc) = i?£(*)iS&(A0 - i££(*)iS&(A0 (60) 

/rl A k' d 4 h" r 

(2^ (2^ (27r)V(A; " k ' ~ k " ] b UiS ^(klriA>^(k'')iS< b (k) 

-YiS< c (k>)Y^(kys? cb (k) 



To check the consistency of the generalised chemical potential ansatz (jHj), we now 
verify that the collision term (I6D1) vanishes provided the leptons and anti-leptons have 
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opposite chemical potentials, fj,~ b = — /i„ 6 . To see this, we first note that fi~ b = —^b — 
—fi a b implies the generalised KMS relation Sf ab {k) = — [e^ k °~^^ Sf cb {k). Using this 
and the fact that the gauge bosons are in thermal equilibrium, which implies A^ (/c) = 

e^A^O), in Eq. §HD yields the analogous relation ff> = - {e /3k °-^) ac ff cb for the 
lepton self energy. This allows us to write 

C?(k)=[iff < {k),e? l »-f , '>] iS<(k) (61) 

Inserting Eq. (j59p leaves the commutator \iS^ l< (k), e^* 0- ^] in flavour space, which is 
easily seen to vanish upon use of the ansatz (1361) together with Eq. (14 ip and n~ b = 
—f^a b - The vanishing of the collision term under the conditions fi~ b = —fi^b means that 
the kinetic equilibrium distribution ()4ip with opposite chemical potentials is indeed a 
stationary solution of the kinetic equation in the limit when only the fast flavour-blind 
gauge interactions are present. 

Furthermore, when the gauge bosons are in equilibrium and the equilibrium deviation 
of the leptons is small and parametrised as in Eq. (|43l) . we can approximate the collision 
term (|60|) as 

Cf ah {k) ^g 2 j ( ^ff 4 (^W* - k ' - k ") ( 62 ) 

YiSS £ab (k')Y [iA^(kyS< b (k)-iAi^(k'')iS> b (k)] 

+ [Yis> aa (k')Y^> Allu (k") -ris< a (k')Y^(k")} m ab (k)} . 

For this expression, we note that the terms in square brackets are odd under a change of 
sign of the momenta, since to leading order in deviations from equilibrium, we may sub- 
stitute the equilibrium distributions and make use of the fact that iA^^(A;) = iA^,(— k) 
and iS^a(k) = ^^(—k). Hence, after performing the k° integration of the collision 
term, the same sign contributions to the equations (J47J) and (|52|) for 5nf and 5nJ occur 
due to a cancellation of the relative sign in Eq. (|47|) . 

However, if we substituted tree- level propagators in Eq. (|62l) . this collision term would 
vanish, since all the three particles involved are massless. It is therefore necessary to 
account for thermal masses and for finite width effects, that relax the zero-temperature 
on-shell conditions. When employing finite width propagators, analytical simplifications 
of the collision integral due to on-shell 5-functions no longer apply. In the present 
work, we therefore do not perform collision integrals that vanish for zero-temperature 
propagators explicitly. Rather, we discuss their general form and give estimates, while 
relegating more precise numerical evaluations to future studies. 

Of particular interest within the collision term ( )62l) are contributions for which 
sign(fc"°) = — sign(fc'°) = sign(/c°). These are allowed when we account for the finite 
width in the spectral functions (for both, £ and the gauge fields A) and they corre- 
spond to lepton anti-lepton pair creation and annihilation processes. After performing 
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the integrations and the Dirac trace, lepton- and antilepton contributions are identical. 
Therefore, we may parametrise the flavour blind contribution to the collision term in 
Eq. (g7D by 

00,0 

± \ tr / ir / || (C + c2) = ~ rbl {6 < b + 6n ^ ' (63) 
0,— 00 

which leads to the corresponding term in Eq. (I52p . Here 5nf ab has been factored out by 
substituting Eqs. flSD and (jS]). Note that by use of Eqs. (S3D and flU} this defines T bl , 
which therefore may readily be evaluated within a more detailed numerical study. For 
the purposes of the present work, we estimate T bl ~ g%T, where gv is the SU(2)l gauge 
coupling and where the additional factor of g\ compared to the tree-level matrix element 



arises from the finite- width effects 22] 



The fact that the flavour-blind collision terms for 5nJ and 5nJ are of the same 
sign also implies that in the absence of additional flavour-sensitive effects, bn[ — 5nJ is 
conserved, as it is required for the ansatz of generalised chemical potentials (l4Tj) to be 
valid. 

3.4 Flavour Sensitive Interactions 

We now turn to the active lepton Yukawa couplings. These contribute to the self-energy 
of the left-handed leptons as 

*£?(*) = hlh db J 0^0^W A S\k -k>- k") iSH d (k')iA^(k") . (64) 
To linear order in deviations from equilibrium, the collision term is 
CL(k) = ^: c (k)iS< b (k) - itl(k)iSUk) « / 0^0^( 2 ^ 54 ( k ~ k '- k ") ( 65 ) 
x {hlh de i5S Rcd (k') [iA>(k")iS< b (k) ~ iA<(k")iSf eb (k)] 
+ hlh de [iS> cd (k')iA>(k") - iS< cd (k')iA<(k")} iSS eeb (k) 

;\J\ (2n) A 5 A (k-k'-k") 

x [hlh db \5S Kcd {k') [iA^kyS^^-iA^kyS^ik)} 

+ hlh ce [iS^>(k')iA>(k") - iS^<(k')iA<(k")} iSS eeb (k)} ■ 

Again, the leading thermal corrections to the propagators should be employed, since 
at tree-level, this integral is vanishing for kinematic reasons. We have to distinguish 
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two relevant kinematic situations: First, when sign(A;°) = — sign(/c'°), the collision term 
corresponds to pair creation or annihilation of a left- and a right-handed Standard Model 
lepton. Second, when sign(/c°) = sign(k'°) the left- and right-handed leptons scatter from 
a Higgs boson. Again, both configurations are only possible due to the finite width of 
the spectral functions of £, R and <fi. We summarise both contributions to Eq. ( H7|) by 
writing 

00,0 

r " = ± i tr / w I H ( c - w + <&<*>) (66) 

0,— 00 



=r an [\tfh\aM* + SngjtfhU + hl5nl cd h db + ^MlA 
+r sc ([tfhlML + Snfl[rfhU - hlMtdhdb - h) ad 5n^ dc h 



d> 



For later use we note that for the right handed leptons, we have the corresponding flavour 
sensitive scattering rate 



rZ =r an {[hh%5nt b + 6n*Jhh% + h ac 5nf c A b + haMlA) (67) 

+r sc ([hrflaMicb + $n^ ac [hh% b - h ac 8n% d h\ h - h ad 5nflhl b ) . 



We estimate the factors r an and T sc as ~ g^T, which is again due to the finite width 
of £ and at finite temperature. (The U(l)y contribution is smaller because of the 
smaller gauge coupling and the smaller number of gauge bosons). We recall that within 
these expressions, 5n e and the second index of the coupling h transform under left- 
handed flavour rotations, while Sn R and the first index of h remain without change. By 
a unitary transformation of the right-handed flavour basis, we may choose the matrix 
htf to be diagonal, which is what we assume here. 

This concludes the derivation of the kinetic equation f )52|) for the number densities. 
An analogous equation (without washout and source terms) holds for the right-handed 
Standard Model leptons. 

3.5 Suppression of Flavour Oscillations 

The largest collision term within the kinetic equations (152|) is T bl = O^T). Close to 
equilibrium, it imposes the constraint 

K b = -<K" 6 • (68) 

This is expected, since in the flavour-blind limit where h ab — > 0, this condition is mani- 
festly invariant with respect to flavour rotations and it reduces to the assumption that 
the lepton charge density of leptons is the same as the lepton charge density of anti- 
leptons. Therefore, this condition is implicitly employed in Ref. [9] as well as in many 
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other kinetic-theory approaches to various problems. Now, due to the ± in the first 
term on the right hand side of the kinetic equations (1521) . a large T bl effectively inhibits 
flavour oscillations, which would be present in the absence of collisions. To see this in 
more detail, consider the toy system of differential equations 

j t 5g + it) = -iAcu8g + (t) - T[5g + (t) + Sg~{t)] , (69a) 

j t Sg-(t) = +iAu8g-(t) - T[5g-(t) + 8g+(t)} . (69b) 

The relevant parameters for flavoured leptogenesis can be estimated as 

r = T bl ~ g\T , Aco ~ h 2 T T < T , (70) 

where h T denotes the r-lepton Yukawa coupling. Since g\ ^> h 2 ., the solutions are linear 
combinations of two eigenmodes with short r s = l/(r + yjT 2 — Au 2 ) ~ 1/(21") and 
long T\ = l/(r — -\/r 2 — Au 2 ) ~ 2T/Au 2 decay times, respectively. The corresponding 
eigenvectors are given by 



^ s ,i = Sg + + 8g w 5g + ± (l t i— J $9 , (71) 

with 

^i = Woe- t/r ^ (72) 

The short mode Sg s ~ 5g + +5g~ is thus damped to zero very rapidly by pair annihilations, 
implying an effective constraint 

5g + ~-(i- i^\ 5g~ . (73) 

Note that the different sign of Aw terms in Eq. ( 1691) is decisive, since it implies that the 
driving term for oscillations in 

j t (5g + (t) - Sg-(t)) = -iAco (6g + (t) + 5g~(t)) (74) 

is damped away, while in the case of same sign 5g + — 5g~ , could have freely oscillated. 
As explained in Section [231 the opposite sign of the Ao; term in Eq. (1521) is a consequence 
of CP invariance at leading order. 

Within the gradient expansion, the first order correction to Eq. (1681) therefore is of 
order Aoj cff /T bl . Since the source terms for the off-diagonal correlations are already 
of first order in gradients, it is justified to use the zeroth order constraint Eq. (1681) 
within our approximations. The long-lived mode describes the damping of flavour 
coherence in the lepton charge density matrix due to flavour-blind interactions. It 
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is much slower compared to the damping rate due to flavour sensitive interactions, 
Au cS /r bl ~ h^g^T C r" ~ O^^tT since h T <C g\. Therefore, we may neglect the 
damping due to flavour-blind interactions, while we keep the direct damping due to 
flavour sensitive processes. While in the case of leptogenesis, we conclude that because 
of Au <C T, flavour oscillations are overdamped and effectively frozen, we note that for 
Au > T, there are damped flavour oscillations. It is interesting to note that even though 
we assume flavour blind interactions, the off-diagonal flavour coherence functions are 
decaying. Such a behaviour, in particular in the oscillatory regime, has been observed 



numerically in Ref. [18 



We emphasise that the conclusion that the oscillations induced by Au are overdamped 
for T ^> Au does not depend on the choice of the flavour basis. To see this, we extend 
g^ 1 to a vector of functions and consider the system of matrix equations 

j t 5g + (t) = -i[u, 8g + (t)} - T[8g+(t) + Sg-(t)] , (75a) 

j t 5g-(t) = +i[w, 8g-(t)] - T[Sg~(t) + 5g + (t)} . (75b) 

Here, T is proportional to the unit matrix and u = u R + u u , where u hl is proportional 
to the unit matrix and u^ b <C T cc for all a,b,c. It then follows that [u^g^it)] = 
[w fl , 5g ± (t)]. By taking the sum of Eqs. (175]) . we again conclude that 5g + (t) + 6g~(t) ~ 
e~ 2r4 . Consequently, the difference of Eqs. (1751) yields 

j t [5g + (t) - Sg-(t)] = + [5g + (t) - 8g~(t)] x O (^j , (76) 

where the right hand side is estimated as the eigenvalues of a matrix with large diagonal 
and small off-diagonal elements. Alternatively, this can be seen by substituting in the 
right hand side of the difference of Eqs. ( J75|) 



[oo,5g + (t) + 5g-(t)] = (5g + (t) - Sg~(t)) x O {oo 2 /T) , (77) 

where an estimate according to Eq. (1751) is made. This confirms the suppression of the 
effect of <^ fl by T bl in a general flavour-basis. 

We can also generalise this discussion to the case of a time-dependent mass basis. In 
order to model this situation, consider the system 

J t 5 9ab(t) = -iAu ab 8g+ b (t) + E ac 8g+ - Sg+E* - T[8g+(t) + 8g~ b (t)] , (78a) 

jfiM = +i&Uab8g- b {t) + E ac 5g a - 5g^ c E cb - T[5g~ b (t) + fy+(i)] . (78b) 

In the limit V 3> Au ab , we may find an approximate solution by imposing 5g + = —5g~. 
This leads to 

j t ($9+ b (t) ~ fe(*)) = K(*) " 8g*{t), E] , (79) 
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which is solved by the unitary evolution 

59 + ab {t) - 8g- h {t) = (Te~ S4 ) {6g+(t = 0) - 8g*(t = 0)) (fe Et ) , (80) 

where T implies the time-ordered exponential. Therefore, the freezing of flavour oscil- 
lations also persists when we account for the time dependence of the mass basis. From 
above equation, we recover Eq. fT76|) by undoing the flavour rotation, that is by left 
multiplication by U and right multiplication by W. 

3.6 Kinetic Equations for Left and Right Handed Number Den- 
sities 

We now define the charge number density matrix as 

qtab = 5n+ ab - 5nJ ab . (81) 

Imposing that fast (compared to the interactions accounted for in W a b and S a b) pair 
creating and annihilating interactions enforce the constraint 

we can take the linear combinations from Eq. (1521) that solve for the charge density 
matrix (1811). For the flavour- sensitive interactions, define 



Tl h =r an ([h)h] ac q idb + q\ ac [h)hU - hlq Rcd h db - h[ d ql dc h cb ) (83) 

+r sc ([tfh] ac q ecb + q\ ac [h ] h] cb - hl c q Rcd h db - h\ d q Kdc h cb ) . 



Using this and the results of the previous sections, we obtain the kinetic equations (J2J) 
which we repeat here for completeness: 

—F. = 2_^ [q^ac^ch — ^acqicb ~ Wacqtcb ~ qtacWdb] + ^S a b ~ ^ £ a b ■ (84) 

c 

Note that similar to the toy system of equations the flavour-blind term drops out in the 
equation for qe ab , while it is consistent to neglect the Au^ term, which would multiply 
Snj' ab + 5nJ ab in this equation, which is strongly damped. This holds in an arbitrary, time- 
independent basis in flavour space, where the H terms are absent. In a time-dependent 
basis such as the basis where q^b is diagonal, the S terms are introduced to account for 
the time-dependent basis rotation. 

For the right-handed leptons, there is the analogous equation 



dqRab ^fl 



dr] 
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Iro* (85) 



with 

rL, =r an {[hh^] ac q Rcb + ql ac [hh% - h ac q £cd hl - h ad q\ dc h\ b ) (86) 

+r sc (\hh ] ] ac q Kch + qli ac [hh% b - h ac q £cd h db - h ad q\ dc h\ h 



Similar results have been obtained earlier within an approach that makes use of the 
density matrix in an occupation number basis. In its details, the equation for the differ- 



ence between lepton and anti-lepton densities in Ref. [13| exhibits however differences to 
our kinetic equation (I84|) . It is not clear whether the lepton charge densities in Ref. 13 
should correspond to our qi^ (which is the difference of the lepton density and the 
transpose of the anti-lepton density) or to the difference of the lepton density and the 



anti-lepton density. In the former case, the flavour oscillations frequencies in Ref. [13 
should have opposite signs for particle and transposed antiparticle modes, if they were 
to agree with our result obtained within the CTP formalism. This is apparently not the 



situation within the equation for the lepton charge density in Ref. [13j. In the latter 
case, as it follows from Eq. ( 13 9 p and the discussion at the end of Section 12. 3\ within the 
CTP formalism the washout and source matrices for the lepton and anti-lepton densities 
are transposed (or complex conjugated, as these matrices are Hermitian) with respect 



to each other, which is apparently not the case in Ref. [13j. Furthermore, the same 
conclusions on damping of coherence would result from our equations if the charge den- 
sity matrix were defined as 8nj~ ab — 5n~[ ba . Hence, with either interpretation, there is a 



difference between the occupation number formalism result that is derived in Ref. [13 



and the kinetic equation (j2]) derived within the CTP formalism. The phenomenological 



consequence of this can be seen when comparing the present work with Ref. [15j] , where 



the kinetic equations from Ref. [13| are solved numerically. While in the present work, 



we conclude that flavour oscillations effectivel y fr eeze out due to fast pair creation and 



annihilation processes, the results in Refs. [13|, [15| imply that the flavour oscillations are 



important and in particular faster than the flavour-sensitive damping processes. 

4 Solutions to the Flavoured Kinetic Equations 

We are considering a scenario with two lepton flavours and assume that there is one 
dominant Standard Model Yukawa coupling h T . In the basis where the lepton Yukawa 
coupling matrix is diagonal, the matrix h is therefore simply 



h=\ :: ■ (87) 




Provided the [i and e Yukawa- couplings are negligible ti 2 ie T/a(r]) <C H, the realistic 
case with three lepton flavours can be reduced to the present case by separating out 
a linear combination of lepton flavours, for which no asymmetry is produced. This 
corresponds to an unflavoured approximation for the e and /i flavours, cf. the discussion 
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of the unflavoured limit below. We note that Eq. fl2]) is manifestly invariant under 
flavour rotations induced by U, while Eqs. ( l4"6j) and (P4"T|) are not, because the term that 
describes flavour oscillations is given in the diagonal basis. In our approximation, we 
can drop this term, because we have shown in Section [3] that due to the constraints from 
kinetic equilibrium, the time scale for flavour oscillations is suppressed when compared to 
the time scale of decoherence from flavour-sensitive scatterings. We use this freedom of 
choice of a lepton flavour basis and perform the discussion in this section within the time- 
independent basis of charged lepton flavours, which is more transparent than the time- 
dependent basis of the leptonic quasi-particles, that is determined by the diagonalisation 
of <^ fl . Likewise, we present all numerical results in the basis of charged lepton flavoursjj 
In the charged lepton basis, the flavour-sensitive collision terms read 
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In the fully flavoured limit, which we define by the requirement (r an + T fl )/i^ ^> H (note 
that for the present purposes, "fully flavoured" refers to the situation where h^ and h e 
are still assumed to be out-of-equilibrium) , we see that within this setup, the flavour 
sensitive collision terms enforce 
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Qei2 — Qm — <?Ri2 — <?R2i 



0. 



(89) 



This agrees with the expectation that when Standard Model Yukawa couplings are in 
equilibrium, the lepton asymmetries are projected onto the charged lepton basis. 

The processes £ + ~R -H- <fi* (annihilation) and £ + <f) -H- R (scattering) are kinematically 
forbidden when all the three particles involved are massless. At finite temperature, this 
holds no longer true due to effects that at leading order can be either thought of as 
thermal masses and finite widths or as radiation of gauge bosons. The latter point of 



view is taken in Ref. [25[ to calculate T sc . However, important t-channel diagrams are 



not included there and a calculation of T an is not provided. A systematic calculation of 



these rates may be performed along the lines of Ref. 

scope of the current work. Motivated by the partial result of Ref. 

numerical definiteness the estimate 



24j . which is however beyond the 



25 , we take here for 



pan I -ps 



0.7a w T/a(i]) = 1.75 x 10~ 2 T/a{r]) 



(90) 



which should be accurate up to a factor of order unity. Besides, we take here aw = 1/40 
as the weak coupling constant at the scale of about 10 12 GeV. The precise value depends 
on the particular extension of the Standard Model. 



3 Yet, we have used the requirement that computations in both bases must yield the same results as 
a consistency check on the numerical results. 
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Figure 1: Comparison of the relevant rates R = r^ D , H, /i 2 (r an + T sc ) for Scenarios (A) 
and (B). The colour key is H (solid, black), Tj D (dashed, dark blue), r 2 D (dashed, light 
blue), h 2 T (T an + T sc ) with h T = 3 x 1(T 2 , 7 x 10" 3 , 4 x 1(T 3 , 2 x 10" 3 , 10" 3 (from top 
right to bottom left, dotted, red). 

To obtain numerical solutions, we first solve the kinetic equations for the distribution 
of the right-handed neutrinos i\q. They are given in Ref. |9|, and the generalisation from 
the single flavour to the two-flavour case follows by the straightforward replacement 
|Yi| 2 — > ^2 a |Yi a | 2 . We employ this distribution to calculate the washout and the source 
terms within Eq. ([2]). To be specific, we choose thermal initial conditions for N\. For the 
singlet neutrino masses, we choose M\ = 10 12 GeV and M2 = 10 14 GeV. For the Yukawa 
couplings of the right handed neutrinos, we consider two scenarios 



/l.4x 10" 2 1 x 1(T 2 
" I i x KT 1 KT 1 

/l.4x 1(T 2 3 x 10- 3 
i x KT 1 KT 1 



Scenario (A) , 
Scenario (B) . 



(91) 



We vary the Yukawa coupling h T , since this will directly exhibit the dependence of the 
results on the flavour effects, while of course, for a phenomenological study, it would be 
more pertinent to vary the unknown parameters Y and M 12 . 

In Figure [1] we show the interaction rates r fl = /i 2 (r an + T sc ) for different values 
of h T and compare them to the expansion rate of the universe H and to the inverse 
decay rate for the individual flavours a, r^ D = 2W aa as a function of the ratio of M\ 
to the physical temperature, z = a{j])M\jT . Scenario (A) exhibits moderate to strong 
washout in both flavours, where the dominant contributions to the lepton asymmetry 
are generated between z ~ 3 and the point when the lepton asymmetry freezes out, 
Tjd « H . We expect flavour effects to be negligible when r fl < H « rf D during these 



times before freeze out 26], i.e. for h T significantly smaller than 4 x 10 by inspection 



of Figure [TJ On the other hand a fully flavoured description should be applicable when 
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Figure 2: Results for Scenario (A) with h T = 3 x 10~ 2 , 7 x 1CT 3 , 4 x 1CT 3 , 2 x 
1CT 3 , 10~ 3 , 0, from top left to bottom right. The key is Ym (dark blue, solid), Y i2 2 
(light blue, solid), Re[Yn 2 ] (dark red, dotted), Im[Y^i 2 ] (light red, dashed). The densities 
are evaluated in the flavour eigenbasis, which means that the larger the flavour effects 
are (the larger h T is), the smaller are the off-diagonal densities. 
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Figure 3: Results for Scenario (B) with h T = 3 x IO" 2 , 7 x IO" 3 , 4 x IO" 3 , 2 x 
10~ 3 , 10~ 3 , from top left to bottom right. The key is Y iu (dark blue, solid), 1^2 (light 
blue, solid), Re[Y^i 2 ] (dark red, dotted), Im[Y^ 12 ] (light red, dashed). 
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r fl > r" D during the times when the quantitatively relevant contributions to the lepton 
asymmetry are produced, i.e. for h T significantly larger than 7 x 10~ 3 . The numerical 
solutions to the kinetic equations for Scenario (A) are displayed in Figure [2j Shown are 
the absolute values of the entropy normalised asymmetries 



45 9* 1 

where we use g+ = 106.75. Since Yi is Hermitian, we plot the real and imaginary parts of 
Yin- The results confirm our expectations about the validity of the fully flavoured and 
the unflavoured descriptions of leptogenesis. For h T ^ 2 x 10~ 3 the total lepton asymmetry 
ti\Y(\ = Yeu + Yi22 is almost independent of h T , and the off diagonal densities decay away 
after freeze out (z rs 10) only. (Only the initial stages of this decay are visible in the 
plots for h T = 2 x 10~ 3 , 10~ 2 in Figure H] due to the cut off at z = 20. For h T = 0, the 
off-diagonal densities do not decay.) On the other hand in the fully flavoured regime, for 
h T > 3 x 10 -2 , the off diagonal densities are strongly suppressed before freeze-out. This 
confirms that neglecting the off-diagonal densities, an approximation that is commonly 
used in the fully flavoured regime, is indeed justified in this regime. In the intermediate 
regime, where neither of these approximations is valid, the correct lepton asymmetry is 
obtained by solving the full kinetic equation (J2J). 

Scenario (B) is a situation with strong washout due to Y n and weak washout for Y 12 . 
The numerical solutions are displayed in Figure [3j Since r 2 D is now significantly smaller, 
it takes also smaller values of h T before the unflavoured description may be expected to 
be valid, cf. Figure [TJ In the fully flavoured regime, we observe that one of the lepton 
flavours suffers from strong washout while the other one is only weakly washed out. On 
the other hand, in the basis where the source term is diagonal, the flavours apparently 
mix in such a way that both lepton flavours are strongly washed out when flavour effects 



are turned off. This importance of flavour effects for washout is well known [12H14I] , and 
it can be easily understood when we recall how the fully flavoured and the unflavoured 
regimes are described. 

First, in the fully flavoured case, densities that are off-diagonal in the flavour basis 
undergo fast damping through flavour-sensitive interactions. As a result, there are two 
washout rates that are proportional to jl^l and \Y± 2 \, respectively. Second, in the 
unflavoured regime, it is convenient to bring Y to a triangular form through 

Y* = Y lb V ba , (93) 

where 

V = - l ( Y22 Y2 *A . (94) 

V\Y2i\ 2 + \Y 22 \ 2 \ -Y 21 Y 2 \ J 

In the triangular basis, a lepton asymmetry is only produced for the linear combination 

V | J 21 1 T 1 1 22 1 
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The washout rate for this linear combination is proportional to 



I y i2 I' = iv i 2 1 | V | 2 (l^iH^il 2 + \Y l2 \ 2 \Y 22 \ 2 + 2Re [Y^Y^Y^Y^]) . (96) 
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The change of the effective flavour basis in the transition from the unflavoured to the 
flavoured regime therefore explains the apparent change in the washout rates for the 
individual flavours for Scenario (B), that are visible in Figure [3J 

To obtain an estimate of the extent of the intermediate regime in terms of the pa- 
rameter Mi, we now fix the r -Yukawa coupling to h T = 0.007, close to its physical value, 
and vary M\ — > aM\ instead. To keep the effect of the washout term and the source 
term constant, we also scale Yu — > y/aYn and Y\ 2 — > \JolY\i as well as M 2 — > aM 2 . 
This scaling behaviour can be seen when recasting Eq. ((2]) into the form 

zH — — = - < y^ [qiac^cb ~ ^acqtcb ~ W ac q ecb - q ac W cb ] + 2S ab - T^ ab \ . (97) 

The terms on the right hand side now correspond to the physical instead of conformal 
interaction rates per unit volume. Since at fixed z, T/a(r]) — > aT/a(r]), H — > a 2 H and 
qt — > a 3 qg, both sides scale as a 5 , except for the term 1/a x rf, which scales as a 4 . 
Therefore, all the scale-dependence is isolated within the flavour-dependent damping 
rate. 

We solve the kinetic equations for 10 10 GeV < Mi < 2 x 10 14 GeV. Parametrically this 



brings us from a regime where flavour effects are maximal to the unflavoured regime [26 
For comparison we also calculate the lepton asymmetry over this parameter range us- 
ing first the unflavoured approximation (h T = 0) and then using the fully flavoured 
approximation, where the off-diagonal number densities are set to zero throughout the 
calculation. 

The results are shown in Figure SJ We find that both the fully flavoured approx- 
imation and the unflavoured approach lead to accurate predictions of the total lepton 
asymmetry within their expected ranges of validity. The intermediate regime where the 
full kinetic equation needs to be solved ranges from around 5 x 10 11 GeV — 10 13 GeV for 
Scenario (A) and even further for Scenario (B) where the unflavoured behaviour is only 
recovered for Mi > 10 14 GeV. This is because the condition for the unflavoured descrip- 
tion to be valid, r^ D > /i^(r an + T sc ) for a = 1, 2 is only fulfilled for larger values of Mi 
within Scenario (B). Besides, within the flavoured approximation of Scenario (B), the 
flavour a = 2 is only weakly washed out. Therefore, quantitatively relevant contributions 
to the lepton asymmetry arise at earlier times, where rf D /[/i^(r an + r sc )] is enhanced 
compared to this ratio close to freeze-out. As a consequence, the fully flavoured descrip- 
tion of Scenario (B) requires smaller values for Mi when compared to Scenario (A). Note 
that the absolute limits for the validity of the unflavoured or fully flavoured description 
may vary by up to order one factors due to the uncertainty in the overall prefactor of 
r fl . We identify this also as a probable source of the numerical difference between the 
present work and Ref. 26|, where it was found using T an + T sc »5x 10~ 3 T/a(r]) that 



the unflavoured description is valid already for Mi > 5 x 10 11 GeV. 
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Figure 4: Shown is the total lepton asymmetry tr[Ye] = Ym + Yg22 as a function of 
the righthanded neutrino mass M\, for parameters corresponding to Scenarios (A) and 
(B). The result of the full kinetic equations including flavour effects (solid blue line) is 
compared to the results of the unflavoured approximation h T = (red, dotted) and to the 
results from a fully flavoured approximation that neglects off-diagonal flavour excitations 
(green, dashed). 
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5 Conclusions 

Using the CTP formalism, we have derived and solved kinetic equations that describe 
flavoured leptogenesis. Our results allow for systematic calculations of the lepton asym- 
metry within the intermediate regime, that is neither fully flavoured nor unflavoured, 
and where off-diagonal correlations between the lepton densities are of importance. The 
CTP framework proves particularly suitable for this problem, since off-diagonal densities 
are straightforwardly implemented within the two-point Green functions. 

So far, kinetic equations that describe flavoured leptogenesis within the intermedi- 
ate regime have only been available as extrapolations from a toy model description of 



number density matrices in the occupation number formalism (13|. The importance of 



a more systematic derivation of these equations has been emphasised for example in 



Ref. 26]. Our main result, Eq. fl2]), that is derived within the CTP framework with the 



Lagrangian ([T]) as the starting point, turns out to resemble the corresponding equation 



of Ref. [13[ in many details of its flavour structure, but it also exhibits qualitative differ- 
ences. The main improvements provided with the present work may be summarised as 
follows: 

• We derive the impact of the thermal dispersion relations on the kinetic equations 
for the left-handed leptons. This confirms the expectation that the effects of these 
dispersion relations share some qualitative features with flavour oscillations in- 
duced by tree-level mass terms and allows for quantitative predictions. Our results 
might also be useful for describing the dynamics of neutrino flavours in interacting 
backgrounds. 

• We find that fast pair creation and annihilation processes through gauge inter- 
actions effectively overdamp flavour oscillations. This is a qualitatively distinct 



feature from the results of Refs. 13, 15 



The washout, source and damping terms in the kinetic equations are derived from 
first principles. They correspond to collision terms which we explicitly present 
in the form of integrals. Some of these integrals crucially depend on finite width 
effects, which make their evaluation difficult. For the purpose of our numerical 
examples, we make only estimates for those collision terms that strongly depend 
on finite- width effects. Yet, these collision terms are well-defined, and their quan- 
titatively accurate evaluation may be the subject of future work, possibly using 



the methods of Ref. 24 . 



For a more accurate prediction of the lepton asymmetry, the present analysis has to 
be supplemented by a number of improvements, which vary in how detailed they have 
been discussed in the literature yet and in how straightforwardly the present work can 



be generalised to include them. First, there are the so-called spectator processes [2711281]. 
which are transitions induced by Yukawa couplings and by strong and weak sphalerons 
and which transfer charges between the £ a and <ft to other particles of the Standard 
Model. Note that also the scatterings induced by h, that transfer charges to the charged 
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right-handed leptons belong to this category. Depending on whether the additional 
interactions are fully equilibrated, out-of equilibrium or in an intermediate regime, the 
kinetic equations have either to be supplemented by algebraic constraints for the various 
charges, or the network of equations has to be extended in a way that is similar to 
how we account for the charged right-handed leptons. Conceptually more interesting 
and challenging is the systematic inclusion of thermal effects. In the present work, we 
have noted that the finite width and thermal mass effects allow for certain three-body 
processes that are kinematically forbidden in the vacuum. Again, the CTP formalism 
bears the potential to account for thermal effects more systematically and may hence 



serve to confirm or to extend earlier results on these effects [29|. Therefore, important 
improvements remain to be incorporated in order to calculate the lepton asymmetry 
of the Universe to a good accuracy and with a quantifyable account of the theoretical 
uncertainty. In order to achieve this goal, the systematic computation of the flavour 
effects from first principles as presented in this work may serve as a building block. 
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A Pole-Mass Equation and Finite Width Propaga- 
tors 

In this Appendix, we show how close to equilibrium, the constraint equation (|12a[) and 
the equations for the retarded and advanced propagators filial) can be solved consistently 
to first order in gradients. The first order corrections account for modified dispersion 
relations and for finite widths. 

In thermal equilibrium, the collision term is vanishing on the right-hand side of the 
constraint equation (JT2aJ. Since the collision term is already first order in gradients, we 
may therefore neglect it within the constraint equations when being close to equilibrium. 
Furthermore, the propagators and self-energies are approximately flavour-diagonal, such 
that we can express the constraint equation ( |12ap in the simple form 

{?-#,!£<'>}-{#'>, LSf} = 0. (Al) 

The spectral function Sf- is defined through Eq. (J7J). In order to solve the pole-mass 
equation (II lap , it is useful to introduce 

n v. ~\ 

(A.2) 
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where Yf t = |tr7 M ^. This facilitates the inversion of propagators within the four com- 
ponent formalism in analogy with the doubling of degrees of freedom within the Weyl- 
fermion propagators that is familiar from the in-out framework 30[. We obtain (c/. 22]) 
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Note again that we have used simplifications due to the flavour-diagonal structure in 
equilibrium. Similarly, we find the Hermitian propagator 
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Now, if we use that in equilibrium 



^(fc°)/;r(k) - 0(-*°)(i«* - fz~m = s ab - 



,/3fc0 + J 



and the KMS relation to substitute 



i^ k ° 
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(A.5) 



(A.6) 



we find that Eq. (I36aj) with Sf- given by Eq. (IA.3J) indeed solves the constraint equa- 
tion (lA.lj) in equilibrium. In a similar way, the same observation holds for Eq. ( 136b)) . 
A related discussion can be found in Ref. 19|. Note that relation (1A.6J) establishes the 
connection between the finite width and the collision term ( TT3"|) . that controls how fast a 
small perturbation in the lepton density relaxes to its equilibrium value. The zero-width 
approximation ( 137)) follows from the result ( 1A.3|) in the limit Yif ■ {k — Hf) — > 0. 

Close to the poles, where ($ — ^ aa ) 2 = 0, the Hermitian propagator Sf is suppressed 
compared to the spectral function Sf by an additional factor Oijl^/k). To first order 
in the gradient expansion and in the narrow width limit, where Y*f ■ {k — T,f) <C k° , we 
can therefore neglect the terms involving Sf in Eqs. ()T2"j) . 

Plugging the ansatz (12 ip into the constraint equation ( lA.lj) and neglecting the term 
{ff J> ASe 1 } leads us to 

gtfk* + h\k\g<f - \ {<T bl + <T fl , 9 <n - \ {sign(^)(,- bl + f) - ^ - , fl , g<>] = , 

(A.7a) 

g<fk° + h\k\g< > - \ {^ + f, g< 3 >} - \ {sign(^)(<T bl + f) - ^ - , fl , g< >] = . 

(A.7b) 
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When neglecting the hole modes, for leptons (k° > 0) the helicity h — — 1 is negative, 
while for anti-leptons (k° < 0) the helicity h — 1 is positive. In conjunction with the 
constraint ( 122]) . this implies within the flavour-diagonal basis the dispersion relations 



k° = ± 



|k|+? 



i)i 



V 



+ 4) 



(A.8) 



for 57^". Note that in the present case, <^ bl ~ g\, while <^ a ~ [V/i] aa . Since g\ ^> [^/i] aa , 
the expressions for the dispersion relations (IA.8J) . that are accurate to order g%, are not 
reliable to order [h^h] aa in case g\ £ [h^h] aa . However, since the flavour-blind terms are 
universal, the differences between the dispersion relations for different i,j are nonetheless 
accurate to order [h^h] aa . 



B Thermal Lepton Dispersion Relation Induced by 
the Right Handed Neutrino 

In this Appendix, we calculate the thermal correction to the dispersion relation of the left 
handed lepton £ induced by its Yukawa coupling Y to the heavy right handed neutrino 
Ni and the Higgs boson <fi. The leading order contribution comes from the one-loop 
wave-function correction where N\ and <fi are running in the loop. A similar calculation 
has been performed in Ref. 23| for a massless fermion in the loop and in Ref. 3l| for 
the light massive case. In the latter calculation, instead of the Higgs field, there is a 
massless gauge boson in the loop, but the intermediate results before taking the limit of 
small fermion mass can be straightforwardly applied to the present situation. 



B.l Decomposition of the Self Energy 

The structure of the self energy and the dressed propagator is restricted by Lorentz 
invariance, which allows us to parametrise both by two invariant functions a and b. 
These can be determined by an explicit calculation of the self energy as performed in 
detail in Ref. 31j . and they can be used to find the poles of the propagator and thereby 
the dispersion relation. We define a and b by: 



h ~2 
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Pi 
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Using Eqs. (I14lll5|16p and (IB.lj) in the rest frame of the thermal bath with u 
we can identify: 

^ = -k°a - b , 

^ = -(|£; |-|k|)a-sign(£;>. 

The poles of the propagators ( 1A.3J) and ( 1A.4I) are given by the equation 



(B.l) 
1,0,0,0), 

(B.2a) 
(B.2b) 

(B.3) 
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which determines the flavoured dispersion relation of the lepton £. Alternatively, the 
dispersion relation is given by the solution of the constraint equations (1A.7|) leading 
to Eq. (IA.8P to first order in gradients (in particular, when assuming a< 1, which is 
justified by perturbatively small coupling constants and loop suppression factors). Sub- 
stituting the relations (IB.2|) we then obtain the following first order dispersion relation 
in the flavour- diagonal basis: 

k° = ±\k\-^(b Baa + b Dbb ) J (B.4) 

where b-Q = U'bU. In accordance with Eqs. ( TT4j 1511 16ft . b can be decomposed into contri- 
butions from flavour blind gauge interaction and flavour sensitive (h- and Y-) Yukawa 
interactions: 

b ab = S ab b hl + hlh cb b^ h + J2 VaYi^V , (B.5) 



with an analogous decomposition for a. We see that according to Eqs. (IB.4J) . ( TB.5J) dif- 



ferent types of contributions to the full dispersion relation are simply additive. There is 
a caveat, however, if the contributions have a nontrivial hierarchy. For instance, in the 
heavy massive case witffl M t > T, we find that &f ,y ~ (T/Nhf while a^' Y ~ (T/Mi) 2 , 
cf. Eq. ( IB. 101) below. Then the term mixing gauge- and Y-induced contributions, 
b bl a{ ~ T 4 /(k 2 M 2 ), would be of leading order in the Y-induced dispersion rela- 
tion ( 1B.4|) instead of b{ , if k 2 /M 2 < 12 g 'j , where C(R) is the quadratic Casimir 



of the fermion representation (see Ref. [23|) and g is the gauge coupling. Using the 
values C(R) = 3/4 and «w = 5 ,2 /(4vr) = 1/40 and assuming |k| ~ T, this condition 
reduces to M\jT > 10. The momentum region |k| ~ T is of particular relevance for 
leptogenesis, since this is where most of the leptons t are present. For z = Mi/T > 10, 
the process of leptogenesis has typically already completed. Thus we will not consider 
the case when these mixed contribution dominate in this Appendix. Instead, we use the 
additive dispersion relation (1B.4J) and consider only the Y-induced contributions in what 
follows. We therefore subsequently also drop the superscript "fl, Y" . 

The one-loop contribution to the thermal lepton self-energy induced by the couplings 
Y is given by 

?S>) = -»M» / J^iA?(p)ft^i(p + *) , (B.6) 

where iAT(p) and iSj^^p+k) are the time-ordered thermal propagators for the Higgs field 
and the heavy right handed neutrino, respectively (cf. Ref. J9(). We have restricted the 
sum over the right-handed neutrino flavours to % — 1, since in the hierarchical case onl 
the lightest right handed neutrino leads to a relevant contribution during leptogenesis 



my 
is@ 



In order to simplify the notation, in this Appendix, we denote by T the physical temperature, in 
contrast to the comoving temperature throughout the remainder of this paper. 

More precisely, this restriction is justified below where we find that in the heavy massive case the 
leading order correction to the dispersion relation is proportional to (T/Mi) 4 . 
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In Ref. 3l|, only the light massive case is discussed where T 3> M\. For leptogenesis 
however, we also need a calculation that is valid at late times, when T <C Mi as well as in 
the intermediate regime which can be treated only numerically. Since our calculation is 
very close to the one in Ref. |3llj we do not repeat the first stepsjj and continue evaluating 
the integrals given in Ref. [31] for the case of a heavy massive fermion in the loop. 

B.2 Heavy Massive Case 

In the heavy massive case T <C M\, we perform an expansion of the logarithmic functions 
L± (|p|) defined in Ref. 3l( and of the neutrino number density np(E) = l/(e E ' T + l) in 



powers of T ' jM\. The energy of the heavy neutrino, E(p) = a/p 2 + Mf, is of order M l5 
and thus the distribution function np{E) is exponentially suppressed for T <C Mi. For 
this reason the fermionic contributions involving L±(|p|) need not be considered. Within 
the functions L±, we count p°, |p| and k°, |k| as order T, such that the expansion up to 
order (T/Mi) 2 yields^ 

T r/i in 8|p| Ikl ( ^ k 2 

Lf(|p|) = 0. (B.7) 



Inserting Eqs. (IB.7P into the expressions for the self-energy (in accordance with Ref. [3l|). 
the remaining integrations reduce to basic integrals and we obtain: 

fr(#£T) = -nX>^ , (B.8a) 

k°T 2 
tr(^f) = -YilY^ ■ (B-8b) 

It is easy to relate these expressions to the F-induced contributions to the Lorentz 
invariant functions a and b decomposed as Eq. (IB. 51) to get 



a i = TTTI72 » ( R9a ) 



Y2M 2 
6i = 0. (B.9b) 



Up to order (T /Mi) 2 the dispersion relation (1B.4I) then simplifies to the massless disper- 
sion, k° = |k|. The calculation up to order (T/M1) 4 requires a more tedious expansion. 



^Thc difference in the calculations affects only the prefactors but not the structure of the integrations. 

rH,Y 

hab 



7 L^(|p|) has to be expanded one order further than L^(|p|) to extract the correct result for tr($K*,' 6 ) 



up to a given order. 
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For the functions a and b we find 

T 2 / 15fc 2 - 4tt 2 
12Mf V 1 + 15JV^ 



«i = 7^777 ( 1 + 77772 ) > (B-lOa) 



47T 2 fc°T 4 

b ^^5Mf (R10b) 

At this order, the non-zero b implies the following dispersion relation for g^ in the 
flavour-diagonal basis: 

where Fj, = | \Ul c Y{ c Y ld U da + Ul c Y{ c Y ld U db ) , and + applies to negative, - to positive 

helicity h. Furthermore, we note that to the leading order (T/Mi) 2 , Eqs. (1B.2J) and (1B.9J) 
imply that 

fl , y 47r 2 |fc°|T 4 _ fl>y fc°T 2 

*>i = 77771 and ^l = "77772- (B.12) 



Note also that a is even in k° and b is odd, such that q' and q' have the correct 
symmetry properties fl35|) . We emphasize that these results are valid forLj |k| < T <C M\. 
In the limit of |k| ^> M\ ^> T, we find another analytical expansion with 

tr(#£ y ) = -n* a ^ , (B.13) 

while the contribution from tr(^^r' 6 ) to the dispersion relation is vanishing in the leading 
order. The dispersion relation is then given by 

fc0 = ± l k l( 1 + ^^)' ( RM ) 

from which an effective thermal mass for leptons £ with large momentum may be ex- 
tracted: 



\Y AT 

y \ I ab\ 1 

2V3 



™heavy = -ITTtT- ( B - 15 ) 



B.3 Light Massive Case 



We consider next the light massive case Mi <C T. Assuming |k| > T and further that 
the thermal corrections to the dispersion relation are small, ||A; | — |k|| <C T, we find 
that to the leading order, b is given by 

bl = IL (J^i n (!^±M) - A (B.16) 

1 16k 2 \2k \k°- k J / y J 



Within this Appendix, we imply by A < B that A is either of order of B or much smaller. 
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Then, by using the assumption of small thermal corrections (in particular a < 1) we 
can use Eq. (IB. 41) to find the dispersion relation 

fc ° = ± l k l( 1 + ^l£)' ( R1? ) 

where we have neglected a logarithmic correction proportional to f Y^-j^p ) m ( ^absik 2 " 

Now we can extract an effective thermal mass for |k| > T by using the relation mX 
k° 2 - k 2 to find 



Y \YJT 

ilight " 2V^ 



<ht = L ^- (B.18) 



For leptogenesis, where most of the leptons are within the momentum region |k| ~ T, 
the results (IB. lip and (1B.17J) imply that the matrix U which diagonalises the lepton 
mass will be approximately constant for T ^> M\ as all contributions to the dispersion 
relations of £ are proportional to T 2 . Then it is time dependent around T « Mi and 
again approximately time independent when T -C Mi, when the contributions induced 
by the couplings Y are suppressed by (T/Mi) 4 . 

B.4 Numerical solutions 

In this Section we solve the F-induced dispersion relation numerically and compare it 
to the analytical approximations ( IB.llj) . (1B.14|) and (IB.17J) . For comparison, we also 
derive the dispersion relation and the thermal mass for the light massive case using hard 
thermal loop (HTL) approximation when k°, |k| ^C T. In this Section we consider only 
the single flavour case and set Y 2 b = \Yi\ 2 = 1 in all plots. Note also that we do not take 
account of relevant contributions from gauge couplings, such that the numerical results 
presented in this Section should be considered as checks of the analytic limits and as a 
study of a toy system in the absence of gauge interactions. 

For the purpose of comparison to the numerical result, we give the analytical ex- 
pression for the dispersion relation for the light massive case Mi <^C T within the HTL 



approximation k , |k| <^C T. Following Ref. [31[ we find 



which leads to 



Tr(#gT ) = VZxJj , (B.19a) 

T 2 fk° 4- lkl\ 
Tr^ST) = YC.Yu.^ln I -^U \ , ( B .19b) 



T 2 ( k» , //t +|k. 

'""leFl 1 -^' 11 !^ 11 ' (&20a) 

6 > - m (m* (^m) - k ° ) ■ (B20b) 
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Figure 5: Shown is the single flavour dispersion relation with k° — |k| as a function of 
|k| in the light massive case Mi <C T. The solid black line corresponds to the numerical 
solution, while the dashed red line to the approximate dispersion relation (IB.21J) and the 
dotted red line to Eq. ( lB~T7j) . 



From this expansion, we see that the suppression of a due to small coupling constants 
and loop suppression factors breaks down in the limit |k| — > 0. We notice that b is given 
by the same expression (IB.16J) as in the light massive case for |k| > T. This is not the 
case for a, however. To extract the thermal mass we expand a and b to the lowest order 
in |k|/|A; | to get a\ = —7^75? and b\ = — ^o- Now we cannot use the dispersion relation 
(IB.4p . since it has been derived assuming that a< 1, which is not the case when |k| — > 0. 
(Using Eq. (IB.22P below, we find that a = 1/3 for |k| = 0). In the single flavour case, 
we find an identical dispersion relation to the one in Ref. 31] : (k° =1= |k|)(l + a) + b = 0, 
which now gives 



fc° — |Yi| 



T 2 



IQk 



±|k| [l-\Yi\ 



r- 



48k ' 



(B.21) 



The thermal mass in the limit |k| — > is then given by 

\Yi\T 



in 



Y 
HTL 



(B.22) 



This is a well-known result within the HTL approximation, see Ref. 23 



In Figure El we compare the approximate dispersion relations (IB. 171) and (IB. 211) for 



the light massive case Mi < T to the exact numerical solution. We confirm that the 
HTL approximation agrees with the exact numerical solution only for small |k|, while 
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Figure 6: Shown is the dispersion relation for the lepton £ with \k° — |k| |/|k| as a function 
of M 1 /T for |k|/T = 0.7, 1.0, 1.3 from dark red (top at left) to light red (bottom at left). 
The dashed line correspond to the analytic approximation (IB. lip . 
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Figure 7: Shown is the dispersion relation with k° — |k| as a function of |k| in the heavy 
massive case with Mi/T = 8. The solid black line corresponds to the numerical solution, 
while the red dashed and dotted lines are the approximate dispersion relations (IB. lip 
and ( lB.14p . Also shown is the light cone k° = |k| (horizontal line). 
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Figure 8: Shown numerically obtained dispersion relations with k° — |k| as a function 
of |k|, for Mi/T = 0.7,0.9,1.0,1.2,1.5,2.0 from dark red (top) to light red (bottom). 
Also shown is the case with M\/T = (black, top) and the light cone k° = |k| (black, 
middle). 



approximating the dispersion by the thermal mass (1B.18J) is reasonably accurate in the 
region |k| > T. 

We next turn to heavy massive case M\ ^> T. In Figure [6] we compare the large 
M\jT limit of the approximate analytical dispersion relation (IB. Ill) to the numerical 
solutions to find a reasonably good agreement for M\/T > 10. In addition, these plots 
clearly confirm the asymptotic (Mi/T) 4 behaviour. 

Moreover, we notice that the approximate dispersion relations ( IB. Ill) and ( 1B.14J) give 
|A; | < |k| for small |k| and \k°\ > |k| for large |k|, so we conclude that in between their 
regions of validity the dispersion curve must cross the light cone k° = |k|, implying an 
intermediate region of superluminal group velocity v g = 4^4 > 1, which might indicate 
dissipative effects in this momentum region. This interpolation is indeed confirmed by the 
exact numerical solution of the dispersion relation. In Figure [7] we show the approximate 
dispersion relations together with exact numerical solution for M\jT = 8. We see that 
the analytic approximations are fairly accurate in the regions |k| < T and |k| ^> T and 
the crossing of the light cone takes place at the momentum |k| ~ 10T. 

The region in between large and small Mi/T can be treated only numerically. In 
Figure |8] we show the numerical solution for the dispersion relation for several values of 
Mi/T. We observe that for M\/T > 0.9 the dispersion curves cross the light cone at 
larger |k| for larger M\. It would be of interest to investigate how finite-width effects 
and thermal masses for the neutrino and the Higgs fields affect the present analysis. 



40 



References 

[1] E. W. Kolb and S. Wolfram, "Baryon Number Generation In The Early Universe," 
Nucl. Phys. B 172 (1980) 224 [Erratum-ibid. B 195 (1982) 542]. 

[2] J. S. Schwinger, "Brownian motion of a quantum oscillator," J. Math. Phys. 2 
(1961) 407. 

[3] L. V. Keldysh, "Diagram technique for nonequilibrium processes," Zh. Eksp. Teor. 
Fiz. 47 (1964) 1515 [Sov. Phys. JETP 20 (1965) 1018]. 

[4] W. Buchmuller and S. Fredenhagen, "Quantum mechanics of baryogenesis," Phys. 
Lett. B 483 (2000) 217 |arXiv:hep-ph/0004145| . 



[5] A. De Simone and A. Riotto, "Quantum Boltzmann Equations and Leptogenesis," 
JCAP 0708 (2007) 002 |arXiv:hep-ph/0703175] . 

[6] M. Garny, A. Hohenegger, A. Kartavtsev and M. Lindner, "Systematic approach 
to leptogenesis in nonequilibrium QFT: vertex contribution to the CP-violating 
parameter," Phys. Rev. D 80 (2009) 125027 jarXiv:0909. 15591 [hep-ph]]. 

[7] M. Garny, A. Hohenegger, A. Kartavtsev and M. Lindner, "Systematic approach to 
leptogenesis in nonequilibrium QFT: self-energy contribution to the CP-violating 
parameter," Phys. Rev. D 81, 085027 (2010) |arXiv:0911.4T22l [hep-ph]]. 

A. Anisimov, W. Buchmuller, M. Drewes and S. Mendizabal, "Leptogenesis from 
Quantum Interference in a Thermal Bath," Phys. Rev. Lett. 104, 121102 (2010) 
jarXiv: 1001 .38561 [hep-ph]]. 

[9] M. Beneke, B. Garbrecht, M. Herranen and P. Schwaller, "Finite Number Density 
Corrections to Leptogenesis," larXiv:1002.13 26 [hep-ph]. 

[10] M. Garny, A. Hohenegger and A. Kartavtsev, "Quantum corrections to leptogenesis 
from the gradient expansion," larXiv:l0 05.5385 [hep-ph]. 

[11] T. Endoh, T. Morozumi and Z. h. Xiong, "Primordial lepton family asymmetries in 
seesaw model," Prog. Theor. Phys. Ill (2004) 123 |arXiv:hep-ph/0308276] . 

[12] A. Pilaftsis and T. E. J. Underwood, "Electroweak-scale resonant leptogenesis," 



Phys. Rev. D 72 (2005) 113001 |arXiv:hep-ph7 0506107|. 



[13] A. Abada, S. Davidson, F. X. Josse-Michaux, M. Losada and A. Riotto, "Flavour 
Issues in Leptogenesis," JCAP 0604, 004 (2006) |arXiv:hep-ph/0601083] . 

[14] E. Nardi, Y. Nir, E. Roulet and J. Racker, "The importance of flavor in leptogene- 
sis," JHEP 0601, 164 (2006) |arXiv:hep-ph/0601084] . 



41 



[15] A. De Simone and A. Riotto, "On the impact of flavour oscillations in leptogenesis," 
JCAP 0702 (2007) 005 |arXiv:hep-ph/0611357] . 

[16] T. Konstandin, T. Prokopec and M. G. Schmidt, "Kinetic description of fermion 
flavor mixing and CP-violating sources for Nucl. Phys. B 716, 373 (2005) 
|arXiv:hep-ph/0410135] . 



[17] T. Konstandin, T. Prokopec, M. G. Schmidt and M. Seco, "MSSM electroweak 
baryogenesis and flavour mixing in transport equations," Nucl. Phys. B 738 (2006) 
1 [arXiv:hep-ph/0505103] . 



[18] V. Cirigliano, C. Lee, M. J. Ramsey-M usolf and S. Tulin, "Flavored Quantum Boltz- 
mann Equations," Phys. Rev. D 81 (2010) 103503 |arXiv:0912.3523l [hep-ph]]. 

[19] T. Prokopec, M. G. Schmidt and S. Weinstock, "Transport equations for chiral 
fermions to order h-bar and electroweak baryogenesis," Annals Phys. 314 (2004) 
208 |arXiv:hep-ph/0312110] . 

[20] T. Prokopec, M. G. Schmidt and S. Weinstock, "Transport equations for chiral 
fermions to order h-bar and electroweak baryogenesis. II," Annals Phys. 314 (2004) 
267 [arXiv:hep-ph/0406140] . 



[21] E. Calzetta and B. L. Hu, "Nonequilibrium Quantum Fields: Closed Time Path 
Effective Action, Wigner Function and Boltzmann Equation," Phys. Rev. D 37 

(1988) 2878. 

[22] B. Garbrecht and T. Konstandin, "Separation of Equilibration Time-Scales in the 
Gradient Expansion," Phys. Rev. D 79 (2009) 085003 |arXiv:0810.4"0T6l [hep-ph]]. 

[23] H. A. Weldon, "Effective Fermion Masses Of Order Gt In High Temperature Gauge 
Theories With Exact Chiral Invariance," Phys. Rev. D 26 (1982) 2789. 

[24] P. B. Arnold, G. D. Moore and L. G. Yaffe, "Transport coefficients in high 
temperature gauge theories: (I) Leading- log results," JHEP 0011 (2000) 001 
|arXiv:hep-ph/0010177| . 

[25] M. Joyce, T. Prokopec and N. Turok, "Nonlocal electroweak baryogenesis. Part 1: 
Thin wall regime," Phys. Rev. D 53 (1996) 2930 |arXiv:hep-ph/9410281| . 

[26] S. Blanchet, P. Di Bari and G. G. Raffelt, "Quantum Zeno effect and the impact of 
flavor in leptogenesis," JCAP 0703 (2007) 012 [arXiv:hep-ph/0611337] . 

[27] W. Buchmuller and M. Pliimacher, "Spectator processes and baryogenesis," Phys. 



Lett. B 511 (2001) 74 arXiv:hep-ph/0104189|. 



[28] E. Nardi, Y. Nir, J. Racker and E. Roulet, "On Higgs and sphaleron effects during 
the leptogenesis era," JHEP 0601, 068 (2006) |arXiv:hep-ph/05 12052] . 



42 



[29] G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, "Towards a complete 
the ory of thermal leptogen esis in the SM and MSSM," Nucl. Phys. B 685 (2004) 
89 |arXiv:hep-ph/0310123] . 

[30] A. Denner, H. Eck, O. Hahn and J. Kublbeck, "Compact Feynman rules for Majo- 
rana fermions," Phys. Lett. B 291 (1992) 278. 

[31] E. Petitgirard, "Massive fermion dispersion relation at finite temperature," Z. Phys. 
C 54 (1992) 673. 



43 



